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Abstract. Search for specific chemical signatures of intermittent dissipation of turbulence in diffuse molecular 
clouds. We observed HCO + (1-0) lines and the two lowest rotational transitions of 12 CO and 13 CO with an excep- 
tional signal-to-noise ratio in the translucent environment of low mass dense cores, where turbulence dissipation 
■ is expected to take place. Some of the observed positions belong to a new kind of small scale structures identified 

f — ' in CO (1-0) maps of these environments as the locus of non-Gaussian velocity shears in the statistics of their 

turbulent velocity field i.e. singular regions generated by the intermittent dissipation of turbulence. We report 
the detection of broad HCO + (1-0) lines (lOmK < T£ < 0.5 K). The interpretation of 10 of the HCO+ velocity 
components is conducted in conjunction with that of the associated optically thin 13 CO emission. The derived 
^—^ ■ HCO + column densities span a broad range, 10 11 < iV(HCO + )/Au < 4 x 10 12 cm- 2 /km s _1 , and the inferred 



HCO+ abundances, 2 x 10~ IU < X(HCO+) < HT 8 , are more than one order of magnitude above those produced 
by steady-state chemistry in gas weakly shielded from UV photons, even at large densities. We compare our 
results with the predictions of non-equilibrium chemistry, swiftly triggered in bursts of turbulence dissipation and 
Oh' followed by a slow thermal and chemical relaxation phase, assumed isobaric. The set of values derived from the 

observations, i.e. large HCO + abundances, temperatures in the range of 100-200 K and densities in the range 
100-T0 3 cm -3 , unambiguously belongs to the relaxation phase. The kinematic properties of the gas suggest in 
turn that the observed HCO + line emission results from a space-time average in the beam of the whole cycle 
followed by the gas and that the chemical enrichment is made at the expense of the non-thermal energy. Last, 
we show that the "warm chemistry" signature (i.e. large abundances of HCO + , CH + , H2O and OH) acquired by 
the gas within a few hundred years, the duration of the impulsive chemical enrichment, is kept over more than 
thousand years. During the relaxation phase, the H2O/OH abundance ratio stays close to the value measured in 
' diffuse gas by the SWAS satellite, while the OH/HCO + ratio increases by more than one order of magnitude. 
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1. Introduction 

Parsec scale maps of molecular clouds with only moderate 
star formation provide a comprehensive view of the envi- 
ronment of dense cores before their disruption by young 
star activity. Whether the tracer is extinction, molecu- 
lar lines or dust thermal emission, large scale maps reveal 
long and massive filaments of gas and dust. Dense cores are 
most often nested in these filaments, which are likely their 
site of formation. The open questions raised by star for- 
mation thus include those on dense core and filament for- 



mation. Filaments being by essence structures at least one 
order of magnitude longer than thick, their observational 
study requires maps with large dynamic range. These ob- 
servations, now available, open new perspectives: filaments 
are ubiquitous and span a broad range of mass per unit 



length (e.g. Abergel et al. (1994) Mizuno et al. (1995) 



and Haikala et al. (2004) for massive filaments and 



Hily-Blant (2004)[ |Falgarone, Phillips fe Walker (1991)1 
Falgarone, Pety fc Phillips (2001)|for most tenuous ones). 
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Recent observations of the environment of low mass 
dense cores in 12 CO and 13 CO rotational transitions, have 
also brought to light previously hidden velocity struc- 
tures. The turbulent velocity field of core environments, 
alike turbulent gas, exhibits more large velocity shears 
at small scale than anticipated from a Gaussian distri- 
bution and the spatial distribution of the non-Gaussian 
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occurrences of velocity shears delineates a new kind of 
small scale structure in these core environments: their 
locus forms a network of narrow filaments, not coincid- 
ing with those of dense gas <|Pety fc Falgarone (2003) 



Hily-Blant (2004)| |Hily-Blant, Falgarone fc Pety (20067, 
hereafter Paper IV). It has been shown by Lis et al. (1996) 
that they trace the extrema of the vorticity pro- 
jected in the plane of the sky. As discussed in 
Pety fc Falgarone (2003)[ these velocity structures of large 
vorticity likely trace either fossil shocks which have gen- 
erated long-lived vorticity, or genuine coherent 1 vortices, 
both manifestations of the intcrmittency of turbulence dis- 
sipation, ( |Pety fc Falgarone 2000| ). 

Interestingly, these structures might be the long- 
searched sites of turbulence dissipation in molecular 
clouds 



as proposed by 
after JFPF), and thus the 
chemistry put forward to 



Joulain et al. (1998) (here- 
sites of the elusive warm 
explain the large abun- 



dances of CH + (e.g. |Crane, Lambert fc Scheffer (1995) 



Gredcl, Pineau des Forets fc Federman (2002)[ 



Gredel (1997)^ 



HCO+ ( |Lucas fc Liszt (1996) 

Lucas fc Liszt (2000)|| and now H 2 

( Neufeld et al. (2002)[ |Plume et al. (2004)| ) observed 



in diffuse molecular clouds. Since all the above observa- 
tions are absorption measurements, in the submillimeter, 
radio and visible domains, they suffer from line of sight 
averaging restricted to specific directions and from a lack 
of spatial information regarding the actual drivers of this 
" warm chemistry" . 

The present work is the first attempt at detecting the 
chemical signatures of locations identified in nearby clouds 
as the possible sites of the "warm chemistry" in diffuse 
molecular clouds, i.e. the locus of enhanced dissipation of 
supersonic turbulence. It is part of a series of papers ded- 
icated to the properties of this new kind of filaments in 
dense core environments, in the broad perspective of un- 
derstanding how non-thermal turbulent energy is eventu- 
ally converted and lost as dense filaments and cores form, 
and the timescales and spatial scales involved in the non- 
linear thermal and chemical evolution of the dissipative 
structures of turbulence. 

Unfortunately, dissipative structures amount to very 
little mass i.e. the impulsive releases of energy affect only 
a small fraction of the gas at any time and they are dif- 
ficult to characterize because their emission is weak in 
most tracers. A multi-transition analysis of CO observa- 
tions, up to the J = 4 — 3 transition, provides for this 
gas a full range of temperatures from 20K up to 200K 
and densities ranging from 10 4 cm~ 3 down to 100 cm~ 3 , 
(Falgarone & Phillips 1996J. Recent observations in the 
two lowest rotational lines of 12 CO and 13 CO show that 
the gas there is associated spatially, although not co- 
incident, with gas optically thin in the 12 CO(1-0) line 
( Hily-Blant (2004)| |HiTy-Blant fc Falgarone (2006)1 nere " 
after Paper II). 



The present paper reports observations and analy- 
sis of the HCO + (1-0) emission of a few such struc- 
tures previously identified in nearby molecular clouds. 
Each of them is characterized by a different value of 
the local velocity shear measured in the 12 CO lines 
in order to sample different strengths (or stages) of 
the dissipation process. In a companion paper, we 
present the substructure in one of these filaments as 
detected in a mosaic of 13 fields mapped in the 
12 CO(1-0) line with the Plateau de Bure Interferometer 
(PdBI) ||Falgarone, Pety fc Hily-Blant (2006)] herafter 



Paper III). The target fields and the observations are de- 
scribed in Section 2. We detail the line analysis done to 
derive molecular abundances in Section 3. In Section 4, we 
show that they cannot be reconciled with current models 
of steady-state chemistry and that an alternative model 
involving non-equilibrium heating and chemistry reason- 
ably reproduces the main features of our data, although 
this model is by no means final. The discussion in Section 
5 is devoted to the limitations of our study and the impli- 
cations of the interpretation we propose. 

2. Observations 

2.1. The target fields 

The selected targets belong to the environments of 
the dense cores studied by Falgarone et al. (1998) and 
Pety & Falgarone (2003) These cores were selected be- 
cause of their particularly transparent environment which 
helps isolating the gas directly connected to the cores 
and minimizes confusion with unrelated components. One 
is located in the Polaris Flare, a high latitude cloud 
(Hci thausen, Bertoldi fc Bensch 2002| ), and is embedded 
in an environment where the visual extinction derived 
from star counts at high angular resolution (< 1 ar 
cmin) is 0.6 < Ay < 0.8 mag Cambresy et al. (2001) 



1 so-called because their lifetime is significantly longer than 
their period 



The other, L1512 flLee, Myers fc Tafalla 2001| ), is located 
at the edge of the Taurus-Auriga-Perseus complex in 
a region of low average column density at the par- 
sec scale (Nn 2 of a few 10 20 cm -2 , as deduced from 
the low angular resolution 12 CO(J=1-0) observations of 
Ungerechts & Thaddeus (1987)] L1512 appears to be a 
very young dense core because it is one of the least cen- 
trally condensed among nearby dense cores mapped in 
the submillimeter emission f Ward-Thompson et al. 1994| ) 
and has no signpost of infall motions toward a cen- 
tral object, such as the so-called blue line asymmetry 
( |Gregersen fc Evans 2000] ). It is also one of the denses core 
where the observed molecular linewidths are the narrow- 
est, yet not purely thermal ( |Fuller fc Myers 1993| ). 

Four positions have been observed in the Polaris envi- 
ronment and one in that of L1512. They are shown in 
the 12 CO(J=2-l) channel maps of Figs ^ an d El from 
the IRAM-30m data of |Pety fc Falgarone (2003)] These 
positions sample the full range of small scale velocity 
shears, following the statistics of the velocity field of 
Pety & Falgarone (2003) The observed quantity used in 
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Fig. 1. 12 CO(J=2-l) channel maps of the Polaris field. The 4 positions observed in HCO+ lie at offsets (-150", -75"), 
(-180", 0), (-218", -38"), (-262", -165") relative to (0,0): l H = 123.68, b n = 24.93. They are indicated with crosses 
(positions 1 through 4 respectively). The dense core, not visible in the 12 CO lines, is located along the eastern edge of 
the map, south of offset (0,0). The velocity intervals in km s _1 are given at the bottom of each panel. 
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Fig. 2. 12 CO(J=2-l) channel maps of the L1512 field. The position observed in HCO + is at offset (-55", 0) relative to 
(0,0): RA=05h00m54.5s, Dec=32 39' (1 950). The dense core is centered on the (0,0) position and is about 2 arcmin 
in diameter ( |Lee, Myers fc Tafalla 2001] Pety & Falgarone (2003) I. The velocity intervals in km s _1 are given at the 
bottom of each panel. 



such statistics is the 12 CO line centroid velocity increment while the largest values (positions #2 and #3 in Polaris) 

(CVI) measured for a given lag. The CVI values measured belong to its non-Gaussian tails. The CVI value in L1512 

for a lag of 0.02 pc are given in Table^f 01, each position, also belong to the non-Gaussian tail of the CVI probabil- 

They correspond to local velocity shears (in the plane of ity distribution, but is close to the limit of its Gaussian 

the sky) ranging from 3 to 20 km s _1 pc -1 , the smallest core, 
value in the Polaris field belonging to the Gaussian core 
of the probability distribution of the CVIs in that field, 
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Table 1. HCO + (J=1-0), 13 CO(J=1-0) and C2H(3/2-l/2) line properties derived from Gaussian decompositions. In 
several cases, velocity components with similar centroid and width appear in the HCO + and 13 CO spectra. They are 
labelled A, B, C, D or E depending on the strength of the lines. These are the five cases for which the analysis can 
be completed. The rms noise level a of the spectra is also given. The value of the centroid velocity increment (CVI) 
measured over a lag of 0.02 pc on the 12 CO(1-0) line at each position is also given. 
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Fig. 3. In each of the 4 panels, corresponding to the 4 positions observed in the Polaris field, three lines are displayed: 
13 CO(J=1-0) (thick line) and a multi-component Gaussian fit (thick red curve), 12 CO (either (J=l-0) or (J=2-l), as 
indicated) multiplied to fit the 13 CO temperature scale, and the HCO + (J=1-0) at the bottom with a multi-component 
Gaussian fit (thick red curve). 
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Fig. 5. Same as Fig. |3 for the position in L1512. Note 
the very different antenna temperatures, with respect to 
the Polaris field, for the HCO + line in particular, due to 
the fact that the dense core is sampled by the beam in 
L1512. The velocity scales for the two fields are also very 
different. 



2.2. Observations and results 

The observations were carried out at the lRAM-30m tele- 
scope and spread over three summer sessions (f996, 1998 
and 1999). A few additional 12 CO and 13 CO(J=2-l) lines 
have been obtained in the framework of another project 
during a winter session in 1999. The weather conditions 
never allowed us to observe at high frequency during 
the summers and the lines observed were HCO + (J=1-0), 
13 CO(J=1-0) and C 2 H(3/2-l/2). Frequency switching by 
7.8 MHz was used to subtract the off line emission. The 
autocorrelator had a resolution of 20 kHz providing a ve- 
locity resolution of 0.08 km s" 1 at 89 GHz. The single 




-10 
Velocity (km/s) 



Fig. 6. C2H spectrum at position 2 in Polaris. The line 
parameters derived from the Gaussian fit are: v = —3.5km 
s" 1 , TX=7mK, cr=3mK and Av =3 km s" 1 ). 



sideband noise temperatures of the receivers greatly im- 
proved in 1999 when a new generation of receivers was 
installed. At 3mm, the SSB receiver temperatures were 
in the range 80-130 K depending on the frequency before 
1999 and 50-65 K in 1999. The system temperatures were 
in the range 200-250 K before 1999 and between 150 and 
250 K in 1999. The sideband rejections were estimated to 
be about 20 dB. At 230GHz, the SSB receiver tempera- 
ture was about 90 K with w 13 dB gain rejection. The 
HCO+(J=1-0), 13 CO(J=1-0) and 12 CO(J=2-l) line pro- 
files are shown in Figs. 0- El for the five positions, as well 
as a unique C2H spectrum. 

In L1512, the line of sight crosses the dense core, so 
the spectra comprise the emission of both the dense core 
(the narrow component between 7.1 and 7.3 km s _1 ) and 
that of its environment extending from 6.4 to 8.2 km s" 1 . 
This velocity interval is assigned to the environment of the 
core on the basis of the line profile properties and of the 



large scale channel maps of Falgarone et al. (1998) and 



Falgarone, Pety & Phillips (2001) In Polaris, the lines of 
sight do not intercept the dense core therefore the ob- 
served emission is from the environment exclusively. 
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Fig. 7. Peak antenna temperature versus half-power 
linewidth of all the Gaussian components identified in the 
HCO+(1-0) line profiles and listed in Table □ 



The exceptional signal-to-noise ratio of the line pro- 
files reveals a large variety of complex non-Gaussian line 
profiles, with several broad and weak components. There 
is a general agreement between the total velocity extent 
of the CO and HCO + lines (except for position # 4 in 
Polaris) but the lineshapes are all extremely different. 
The velocity components visible in the 13 CO spectra are 
most often present in the HCO + spectra at the same po- 
sition. The HCO + and 13 CO lines have thus been decom- 
posed into several Gaussians with the goal of ascribing 
a 13 CO(1-0) counterpart to each of the HCO + compo- 
nents. The decomposition of these profiles into Gaussians 
is not straightforward, and certainly not unique, but it is 
a guidance in the analysis of the spectra. Yet, as shown 
in Tabled not all the components found in the Gaussian 
fitting procedure of the HCO + lines have been found suc- 
cessfully in that of the 13 CO lines and vice versa. In most 
cases, the broad HCO + weak emission has been fitted with 
a single Gaussian, although in a few cases, several nar- 
rower Gaussians would also fit to the profile. 

We have restricted our analysis of the HCO + lines to 
the components which are visible in the 13 CO(1-0) spec- 
tra i.e. those for which it exists a 13 CO counterpart with 
comparable centroid velocity and width. They are labelled 
with letters in Tabled Instead of analysing all the profiles 
separately, we give outlines of the solutions by splitting 
the range of line temperatures into five cases (labelled by 
letters A to E). Each case is identified by its HCO + (l- 
0) temperature within a factor of the order of 3, given in 
Tabled Note that in all cases but E, the HCO+/ 13 CO line 
temperature ratio is w 0.1. The weakest cases are the A's 
while the E case approaches the conditions within dense 
cores (noted " dc" in Table . The characteristics of all 
the components are reported in Table Table shows 
that the HCO + line temperatures range between 0.01 and 
0.5 K, a range similar to that found towards £ Oph by 



ing rough anticorrelation appears between the HCO + peak 
intensity and the width of each component where approx- 
imately TX i (HCO+) cx Aw" 2 (Fig.Ej). It will be discussed 
in Section 4.5. 



3. Column densities and abundances 

3.1. HCO+ 

HCO + line intensities have been interpreted in the 
framework of the Large Velocity Gradient (LVG) 
formalism. The collisional excitation due to elec- 
trons has been included in the code following 
Bhattacharyya, Bhattacharyya fc Narayan (1981)[ with 



ionisation degrees x e = 10 4 and x e 



io- 



The HCO+ 



Liszt & Lucas (1994) and in the direction of several ex- 
tragalactic sources by Lucas fe Liszt (1996)| An interest- 



electron rates calculated by Faure & Tennyson (2001) 
have not been used, since they only calculate the 
rates up to J = 2, while for our calculations 
higher excitation levels have been included. We 
have verified that, at least for higher tempera- 
tures, the rates from Faure fc Tennyson (2001)] and 
Bhattac haryya, Bhattacharyya fc Narayan (1981)| do not 
differ significantly. Collisional excitation of HCO + by 
H2 at high temperature is implemented by using the 
cross-sections computed by |Flower (2000)] The excitation 
temperature, optical depth and emergent intensity in the 
J=l-0 transition of HCO + have been computed for H2 
densities ranging between 30 and 10 5 cm" 3 and kinetic 
temperatures T k = 10 and 250 K. Temperatures as large 
as 250 K have been considered because observations of the 
12 CO(J=3-2) and J=4-3 transitions of a few positions in 
the environment of L1512 reveal 12 CO intensities larger in 
these transitions than anticipated from the low J lines if 
the gas is cold and dense. Diluted (rtn 2 = 200 cm" 3 ) and 
warm (Tk = 250 K) gas is a possible solution (Falgarone, 
Pety & Phillips 2001). An upper limit T u = 200 K on 
the kinetic temperature is provided by the width of the 
HCO + velocity components, whose median value is 0.6 
km s" 1 . At T k = 200 K, the thermal velocity of HCO+ is 
only vth — 0.24 km s , and the thermal contribution to 
the linewidth Avth = 0.57 km s" 1 . 

Although weak, the HCO + lines are not optically thin 
throughout the intensity range 0.01-0.5 K. For densities 
lower than a few 10 3 cm -3 , the optical depth approaches 
unity, as the excitation temperature drops towards the 
Cosmic Background temperature. To break the degener- 
acy of the problem, we use the 13 CO(J=1-0) and J=2-l 
lines to constrain the H2 density and temperature domain 
of the solutions and, as said above, we restrict our analy- 
sis to those components of similar kinematic characteris- 
tics in the HCO+ and 13 CO spectra. The 13 CO lines are 
optically thin (or almost thin, in the case of L1512) so 
that the column density per velocity unit is well deter- 
mined and the two constraints of the 13 CO(J=1-0) an- 
tenna temperature and i?(2-l/l-0) line ratio define a nar- 
row track of solutions in the density-temperature plane, 
which nevertheless spans the whole range from cold/dense 
to warm/diffuse gas. Each [HCO + , 13 CO] pair of compo- 
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Table 2. LVG results for the five template velocity components. 13 CO(1-0) line: range of input values for i?(2-l/l-0) 
(col. 2) and T£(l-0) (col. 3), result for the 13 CO column density per unit velocity (col. 4). HCO+(1-0) line: input 
T£(col. 5), lower and upper values for the HCO + column density met along the curves shown in Fig. [S] (col. 6), and 
corresponding range of 13 CO/HCO + ratios (col. 7). 
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nents (i.e. each letter in Col. 4 of Tabled has then to be 
interpreted with the set of values (tih 2 , Ife) allowed by the 
13 CO lines. This method provides determinations of the 
HCO + column density per unit velocity N(HCO + )/Av 
that do not suffer from the lack of knowledge of the den- 
sity, a critical parameter for a line so difficult to excite. 
The results for each (nn 2 ,Tk) set are displayed in Fig. 
where the labels on the curves indicate the gas kinetic tem- 
perature, given only once on the upper curve. The tracks 
appear on that figure, from H2 densities in the range 1— 
2xl0 3 cm -3 at temperatures below 15 K and to densities 
between 100 and 300 cm~ 3 at temperatures above 100 K. 
In the following discussion, we restrict the domain of 
to a range 15 - 20 K to 200 K, the lower value being the 
temperature of a gas of density 1 - 2xl0 3 cm~ 3 , shielded 
from the solar neighbourhood UV field by 0.5 to 1 mag 
\Le Petit, Roueff fc Herbst 2004| ), the upper limit corre- 
sponding to a fully thermal contribution to HCO + line 
widths of 0.57 km s _1 . 

In spite of the large range of kinetic temperatures con- 
sidered, 20 < T k < 200 K, Fig. |H1 shows that the HCO+ 
column density of each component is determined to within 
a factor of a few (see also Table 0). Their values for the 
ensemble of the five template components span more than 
one order of magnitude. 

The main difficulty in the derivation of the HCO + 
abundances from the above set of column densities, is the 
estimate of the relevant hydrogen column density of each 
component. This is so because each component is defined 
only by its velocity and velocity width and contributes to 
only a small fraction of the total column density. In the 
following expression of the HCO + abundance: 



AT(HCO+) = 



Av 



Av- 



1 



(1) 



N ^ H Av' ^ anc ^ nH are derived from the above LVG anal- 
ysis for each case, Av is known from the data. Here, 
njj = 2?ih 2 is the density of protons in a gas where the 
density of atomic hydrogen has been neglected, according 
to PDR models for gas with UV shieldings above 0.2 mag 
( |Le Petit, Roueff fc Herbst 2004| ). The depth L is the only 
independent, unknown and free parameter in the deriva- 
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Fig. 8. HCO + column densities per unit velocity as a func- 
tion of H2 density, for the five template components (A to 
E) for which the LVG analysis is possible (see text). Six 
values of the kinetic temperature are given along the upper 
curve, the same values being marked by the 6 symbols on 
the other curves. Each («h 2 , Tk) set belongs to the narrow 
track of solutions allowed by the CO data (Sect. 3). 



tion of V(HCO + ). An estimate of the depth L is thus 
needed for each velocity component. The channel maps 
in the 12 CO(J=2-l) lines (Figs. Hand |2J) have been used 
to estimate the spatial extent in projection, L, of each 
of the velocity components. The similarity of the velocity 
coverages of the 12 CO and HCO + lines suggests that at 
the scale of the resolution these species are well mixed in 
space too. The lengthscale L characteristic of the struc- 
tures are all of the order of 0.05 pc in projection. The ac- 
tual lengthscale over which the line forms may be longer 
by a factor of a few, less than 10 on statistical grounds 
(then the filaments would be the edge-on projection of 
sheets) or smaller if substructure exists. In the following, 
we adopt L = 0.1 pc for all the components. Note that the 
estimated total column density of individual velocity com- 
ponents thus ranges from N-r = 7.5 x 10 19 cm~ 2 to ten 
times more, or 0.04 < Ay < 0.4, smaller, as expected, 
than the cloud visual opacities discussed in Sections 2 
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is X(C2H) = 2 x 10 8 for an estimated total column den- 
sity A H « 10 20 cm" 2 from n H « 500c m - 3 (Fig. EJ and 
L = 0.1 pc. 

At this position, the same component has ./V(HCO + ) = 
3 x 10 11 cm" 2 , which makes this component similar to the 



10 

n H (cm 3 ) 

Fig. 9. HCO + abundances as a function of the hydrogen 
nuclei density. The five upper curves (thick solid curves) 
are the tracks derived from the observations displayed 
with the same symbols as in Fig. [SJ The kinetic tempera- 
ture varies along each curve and the symbols are located 
at T k =20, 40, 100 and 200 K from right to left. The 
(riH,Tfc) sets are those inferred from the CO data. Thick 
lower curves (red): Results of steady-state chemistry mod- 
els, computed for two different values of the UV shielding, 
A v = 0.5 and 1 mag. Dashed and dotted curves (blue): 
Tracks of non-equilibrium chemical computations along 
two isobaric cooling sequences (see Section 4.4). The ini- 
tial conditions are the same (nn = 100 cm' 3 ) except for 
the shielding from the ISRF: Ay = 0.2 mag (dashed) and 
Ay — 1 mag (dotted). On each of them, the gas tem- 
perature decreases from top to bottom and is given along 
each track to allow a comparison of the LVG and chemical 
solutions in the [X(HCO + ), tir-, T k ] space. 



and 4.2. A velocity coverage Av — 0.5 km s _1 has been 
adopted for all the components as the median value of the 
well-identified components in Tabled The broadest com- 
ponents are likely collections of narrower sub-components. 
The results are shown in Fig.Elfor the five cases analysed. 
The symbols are the same as in Fig. |H1 and on each curve 
they mark the set of kinetic temperatures, T k =20, 40, 
100 and 200K from right to left. According to the previ- 
ous discussion, these abundances may have been over or 
underestimated by a factor 3 at most. We compare these 
determinations to results of chemical models in Section 4. 



3.2. C 2 H 

The width of the C2H line is large and probably corre- 
sponds to the weakest component seen in the HCO + pro- 
file at the same position (case A). With one transition 
detected only, one can simply estimate the column den- 
sity assuming the excitation temperature and low opac- 
ity. For ^3/2,2-1/2,1 = 4.1 x 10~ 7 s _1 and an excita- 
tion temperature of 10 K, the inferred column density is 
A(C2H) = 2 x 10 12 cm~ 2 . The corresponding abundance, 



weakest lines detected by Lucas fc Liszt (2000) in diffuse 
molecular clouds with the PdBI, with an abundance ratio 
C2H/HCO + « 7, consistent with the non-linear relation- 
ship they find between these two species. 



4. Comparison with chemical models 

4.1. Description of the steady-state model for 
chemistry 

Steady-state HCO + abundances have been computed us- 
ing the Meudon PDR code 2 . Models here consider semi- 
infinite clouds of constant density in plane-parallel geom- 
etry with ionizing photons coming from one side. The 
UV spectral density of the photons is scaled according 
to the Draine radiation field (Drai ne~1978[l . The dust ex- 
tinction curve is the galactic one as parameterized by 
Fitzpatrick fc Massa (1990)] The computation proceeds 
from the exterior and penetrates into the cloud. At each 
step, chemical and thermal equilibrium are reached self- 
consistently. Radiative cooling makes use of the on-the- 
spot approximation which considers that an emitted pho- 
ton is either re-absorbed locally or escapes from the cloud, 
and line emissivities depend on the abundances. Radiative 
transfer into lines of H 2 and CO is computed using the ap- 
proximation of Federman, Glassgold & Kwan (1979) The 
depth into the cloud is measured in magnitudes of visual 
extinction (Ay). Therefore, Ay = mag at the edge of 
the illuminated cloud, increasing when moving deeper into 
it. One is then free to consider some particular points into 
the cloud, corresponding to particular visual extinctions. 
We have computed the abundance of HCO + in several 
models with different hydrogen densities, and plot the re- 
sult at two depths in the cloud, Ay = 0.5 and 1 mag (see 
figure 0). Since the total extinction in the direction of the 
observed positions is smaller than 1 mag (see Section 2.1) 
except for the gas in the L1512 dense core, not under con- 
sideration here, it is reasonable to assume that locally the 
UV shielding is also less than 1 mag, hence the two PDR 
cases shown in Fig. |5J 



4.2. Lack of steady-state chemical solutions 

Fig. [5] shows that none of the HCO + abundances de- 
rived from the observations are consistent with steady- 
state chemistry in clouds moderately shielded from the 
ambient UV field, even at high density. The uncertainties 
in the abundance derivations are large and due mostly 
to the unknown extent of the emitting gas along the line 
of sight L. Note, however, that a depth 10 to 100 times 



2 Code available at 
http : / /aristote . obspm . f r/MIS/ index . html 
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larger (i.e. 1 to 10 pc) would be required to bring the 
curves deduced from the observations down to the do- 
main of steady-state solutions. Depths larger than ~ 3pc 
are ruled out by the average column density in the fields, 
Ay < 1. Structures smaller than 0.1 pc are indeed sug- 
gested by the PdBI observations reported in Paper III. We 
thus conclude that the observed HCO + abundances can- 
not be produced by state-of-the-art models of gas phase 
steady-state chemistry in diffuse molecular clouds. 

4.3. Non-equilibrium chemistry 

An alternative scenario relies on the proposal that small 
regions bearing chemical signatures different from those 
of the bulk of gas are permanently generated in inter- 
stellar matter by bursts of dissipation of turbulence. The 
space-time intermittency of turbulence dissipation is such 
that the heating rate due to the local release of suprather- 
mal energy may exceed those driven by the UV radiation 
field. In low density gas and in the UV field of the Solar 
Neighborhood, it is the case as soon as dissipation is con- 
centrated in less than a few % of the volume (see the 
review of |Falgarone 1999| ). 

Intermittent dissipation locally heats the gas far above 
its average equilibrium temperature, over a timescale suf- 
ficient to trigger a specific chemical network driven by 
endothermic reactions or involving large activation bar- 
riers. In the regions heated by non-thermal energy from 
turbulence, the endothermic barrier (AE/k =4640 K) of 
the C + + H2 — > CH + + H reaction is overcome and it 
occurs at a much higher rate than the slow radiative asso- 
ciat ions, C+ + H -> CH+ + hv and C+ + H 2 -> CH+ + 
hv {Black, Hartquist fc Dal garno T978J. Such a scenario 
was first suggested by the good observed correlation be- 
tween CH+ column densities and rotationally excited H 2 



1 Lambert & Danks (1986)). The CH + ions then initiate 
the hydrogenation chain, CH + + H2 — » CHj + H and 
CH^ + H2 — » CH^ + H, which does not proceed fur- 
ther because CH^" recombines with electrons faster than 
it reacts with H2. In the regions heated by non-thermal 
bursts of energy, CH^ thus becomes the most abundant 
ion. As discussed in the specific model of JFPF, it is even 
more abundant than CH + because its formation rate stays 
larger than its destruction rate. It opens a new route of 
formation for HCO + , 



cm 



o -> Hccr 



Ho. 



(2) 



Large amounts of HCO + are thus expected to be produced 
by this chemistry, even in regions weakly shielded from the 
UV field. There, HCO + is, in part, a daughter molecule 
of CH + , although other routes are also active (see Section 
5). 

Following |Moffatt, Kida fc Ohkitani (1994)| the re- 



gions of intermittent dissipation have been modeled in 
JFPF by Burgers vortices threaded by helical magnetic 
fields, an alternative scenario to magneto- hydrodynamic 
(MHD) shocks in the absence of obvious post-shock com- 
pressed layer. The vortex is fed by large scale motions. 



Its lifetime is therefore the turnover time at of these large 
scales, of the order of ten times (or more) its own pe- 
riod, P. Therefore, the characteristics of the vortex are not 
free parameters but imposed by the ambient turbulence. 
Viscous dissipation is concentrated in the layers of large 
velocity shear at the edge of the vortex. In magnetized 
vortices of the CNM, JFPF show that decoupling occurs 
between the neutrals spiralling in and the ions tightly cou- 
pled to steady helical fields. An additional local heating 
rate, due to the friction between the ions and the neutrals, 
is thus present. The chemistry in the diffuse gas reacts 
swiftly to the sharp increase of gas temperature generated 
by turbulence dissipation. In the model of JFPF, relevant 
for the CNM, the period of the vortex is P = 600 yr, its 
lifetime « 10 4 yr, and the minimum time to build up the 
molecular abundances specific of the warm chemistry is 
twarm ~ 200 yr. 

The sequence of events we consider is therefore: dif- 
fuse gas enters such a vortex (or an ensemble of vortices 
braided together), spiralling inward. It is heated and en- 
riched chemically during at least 200 yr while it crosses 
the layers of largest velocity shear and largest ion-neutral 
drift, the active layers of the vortex. Then it enters the 
central regions of the vortex where the temperature drops 
due to the decrease of dissipation. The gas formerly heated 
and chemically enriched, starts cooling down and con- 
densing self-consistently with chemistry. Its radial velocity 
there has become vanishingly small. Eventually, the vortex 
blows-up, after sa 10 yr. 

The thermal and chemical evolution of the gas in the 
active warm layers of the vortex has been computed in 
JFPF. Here, we follow the time-dependent evolution of the 
density, temperature and chemical abundances of the gas, 
once it has left these active layers. An isobaric evolution is 
assumed. We use another code with the same chemical net- 
work and the same cooling functions. The time-dependent 
evolution is computed for a given and constant shielding 
from the UV field. The gas temperature is derived from 
the radiative cooling rate which depends on the chemical 
abundances, the H2 density and temperature. 

The results are displayed in Fig. for four differ- 
ent initial conditions in density and shielding from the 
UV field. Two initial densities are considered nn = 30 
and 100 cm -3 , representative of the cold neutral medium 
(CNM). The local heating rate due to turbulent dissi- 
pation, r^, characterizes the strength of the burst. Two 
values have been chosen, Yd = 10~ 21 ergcm^ 3 s _1 for 
the densest case and = 10 -23 ergcm~ 3 s _1 for the 
other, so that the initial kinetic temperature is close to 
10 3 K in each case. These values correspond to the av- 
erage heating rate due to turbulence dissipation in the 
CNM, T d = 2 x 10- 25 (n H /30cm- 3 ) ergcm _3 s _1 , divided 
by the estimated volume filling factor of the gas affected 
by the burst, or « 2 x 10~ 2 for the weaker type of burst 
and « 10 -3 for the stronger and therefore rarer type (see 
the discussion in JFPF and |Falgarone 1999| ). Two values 
of the shielding from the ambient ISRF arc also considered 
Ay — 0.2 and 1 mag to bracket the conditions prevailing 
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Fig. 10. Time-dependent evolution of (top) the HCO + 
abundance, (middle) the kinetic temperature and (bot- 
tom) the hydrogen density, along the isobaric cooling of 
a parcel of gas formerly enriched chemically in a burst of 
dissipation of turbulent energy. The initial densities are 
TiH,i =30 (dashed and solid) and 100 cm -3 (dotted and 
dot-dashed), and the shieldings from ISRF are Ay — 0.2 
(solid and dotted) and 1 mag (dashed and dot-dashed) 

in the gas under study. The ionisation degree stays close 
to x e = 10~ 4 along the isobaric evolution at Ay = 0.2 
and varies between x e — 10~ 4 and 10~ 5 in the other case, 
following the chemical evolution of the ions. These upper 
and lower values of those used in the LVG code. 

The initial chemical abundances are those produced by 
the warm chemistry in the active layers of the vortex, as 
computed by JFPF. 

It is remarkable that the signatures of the warm chem- 
istry are kept by the gas for more than 10 3 years after the 
gas has escaped the active layers. The curves of Fig. 
illustrate the marked non-linearity of the chemical evolu- 



tion: (1) large differences are visible after only 10 3 years 
between the tracks followed by two parcels of gas with very 
similar initial conditions, (2) the drop of X(KCO + ) is by 
6 orders of magnitude while the density and temperature 
vary by a factor 50 only. 

4.4. Existence of [X(RCO + ),n H ,T k ] tracks crossing 
the observation domain 

We show below that we may have detected the decay- 
ing signatures of the warm chemistry generated by bursts 
of dissipation of turbulence. The time-dependent evolu- 
tion of an enriched parcel of gas once the heating due 
to dissipation of turbulence has dropped, is a curve in 
the [X(HCO + ), «h, Tk] space which depends on the ini- 
tial conditions and UV shielding (see Fig. 110(1 . It is there- 
fore not straightforward to compare it with observational 
results. As an illustration that such an evolution is con- 
sistent with our data, we have drawn two of the tracks 
of Fig. (those corresponding to the largest initial den- 
sity) in the [X(HCO+)-n H ] plane of Fig. M On each of 
them, the gas temperature decreases from top to bottom 
from about 500 to 50 K. In that plane, the LVG solutions 
and the cooling tracks are almost perpendicular. The only 
temperature range where they are consistent are over the 
range 100-300 K: the low extinction track (dashed curve) 
is consistent with the LVG solutions for components A 
and B and the more shielded evolution (dotted curve) for 
all the others (C, D and E). 

The range of values of the HCO + abundances, density, 
kinetic temperature of the part of the isobaric cooling con- 
sistent with the observations (cases A to E) are given in 
Table El for five different sets of initial conditions, includ- 
ing the strength of the dissipation burst which determines 
the local heating rate, r^. 

We have not explored a larger domain of initial condi- 
tions because our goal here is only to illustrate the compat- 
ibility of such a non-equilibrium chemistry with the data. 
A more detailed analysis requires further constraints to 
treat the many non-linearities of the physics and chem- 
istry involved and explore a relevant part of the vast pa- 
rameter space. Moreover, our model, although complex 
in terms of coupled processes, is still too crude. For in- 
stance, we have neglected the initial dynamical expansion 
of the warm gas because the radiative cooling and chemi- 
cal timescales, up to ~ a few thousand years, are smaller 
than the expansion time of the warm gas under the effect 
of pressure gradient, just after the burst of dissipation. 
The timescale for a warm structure of initial thickness Lq 
to double its size under expansion at the sound velocity 
in the ambient cold neutral medium (cs ~ 0.5 km s _1 ) 
is 2x 10 5 (Lo/0-l pc) yr. Would they be much smaller, the 
timescale would be shorter. The Lorentz force acting on 
the ions somewhat delays the expansion at a rate which 
depends on the gas ionisation degree and geometry of the 
helical magnetic fields. 
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Table 3. Sets of non-equilibrium HCO + abundances consistent with HCO + (J=1-0) line observations. Each solution 
is given as a range of values for the set [X(HCO + ), 7iH,Tfc] corresponding to the range of observed values (cases A to 
E). The parameters of the models are the turbulent heating rate T^, the shielding from the UV field A v and the initial 
density nu,i- 
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Nonetheless, the above comparison brings to light a 
few interesting results: all the solutions provided by the 
non-equilibrium chemistry consistent with the observa- 
tions are warmer than 100 K (Table [3J and the inferred 
HCO + abundances are all larger than 2 x 10 -10 and as 
large as « 10 -8 . They are more than one order of magni- 
tude above the values computed at steady-state at densi- 
ties less than 10 3 cm -3 . It is remarkable that these abun- 
dances b racket the value X(HC O + ) = 2 x 10~ 9 deter- 
mined by Lucas &: Liszt (1996)| and |Liszt fe; Lucas (2 000) 



from their observations of the HCO + (J=1-0) line in ab- 
sorption against extragalactic continuum sources. These 
lines of sight sample ordinary interstellar medium of low 
visual extinction (many of the extragalactic background 
sources are visible QSOs) as opposed to dense cores and 
star forming regions. It is thus plausible that their obser- 
vations sample the kind of molecular material we discuss 
here, their inferred abundances being an average along the 
lines of sight of HCO + -poor and HCO + -rich regions. Last, 
results are weakly sensitive to the strength of the dissipa- 
tion burst, as if the actual initial conditions were rapidly 
forgotten. 



4.5. The HCO + enrichment and the strength of the 
velocity shear 

As seen in Table ^ the five observed positions do 
not sample gas with similar kinematic properties: two 
positions (Polaris # 2 and 3) lie on the locus of 
largest velocity shears in the Polaris field (i.e. populat- 
ing the non-Gaussian tails of the probability distribu- 
tions of CVIs, |Pety & Falgarone (2003)| ). The others cor- 
respond to smaller velocity shears. Against all expecta- 
tion (Tabled]), the largest HCO + abundances (D, E and 
dc components) are found at those positions of smallest 
CVIs, while the set of least chemically enriched compo- 
nents (A, B and C) are found at the positions with the 
largest CVIs. Since the larger the velocity shear, the larger 
the dissipation and heating rates and the gas temperature, 
one expects that the chemical enrichment would increase 
with the velocity shear or the CVIs, but it is not what is 
observed. 



These results suggest a broader time sequence, encom- 
passing that discussed in the previous section, in which the 
increase HCO + abundance is a measure of the elapsed dis- 
sipation: the longer the dissipation has been going on, the 
smaller the velocity shear (i.e. the vortex is fading) and the 
larger the quantity of chemically enriched gas produced. 
Such a sequence is also suggested by the rough anticorre- 
lation displayed in Fig. The HCO + column density of 
each component being close to • Au», since the line is not 
very optically thick, it increases as the velocity dispersion 
decreases, as Aw" 1 . The observed quantity of chemically 
enriched gas is produced at the expense of the non-thermal 
kinetic energy, traced by Auj. 

5. Discussion and perspectives 

5.1. The complex space-time average hidden in the 
observed HCO + emission lines. 

The results of the last section seem to be in contradiction 
with the relaxation scenario described above, in which the 
HCO + abundance decreases with time. Another, proba- 
bly related question, is: Why, in regard of the six orders 
of magnitude spanned by the HCO + abundance in the 
proposed non-equilibrium evolution, the observations in 
emission point to a much narrower range of abundances, 
210~ 10 < V(HCO+) < 10~ 8 ? 

This may be understood in a framework in which, sev- 
eral phases of the evolution are present simultaneously in 
the beam, for individual vortices much thinner than the 
beamsize, here 0.015 pc (or ks 3000 AU) while vortices, 
in the JFPF model, have a radius w 15AU. Observations 
sample a mixture of gas: the most diluted and HCO + -rich, 
still in the warm active layers of the vortices, as long as 
the vortex is alive, gas already out of the warm active lay- 
ers and relaxing, up to the densest, HCO + -poor gas. The 
older the burst, the larger the amount of relaxed gas. The 
signal detected in emission is sensitive to a combination of 
the HCO + column density, or simply the abundance X(t) 
assuming the total gas column density being processed 
is constant with time, the density nn(i) that enters the 
line excitation and a term involving the duration of the 
phase, « t, given the broad range of timescales involved. 
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We have computed the time integral of X(t) nn(t), taken 
as a template for the signal in emission, assuming that we 
observe all the phases in the gas accumulated during at 
least ~ 10 4 yr, the lifetime of the vortex. For all the cases, 
this quantity increases linearly with time up to a maxi- 
mum at about 10 3 yr and does not increase further. The 
observations in emission are thus likely to be most sen- 
sitive to the the relaxation phase because it lasts longer 
and involves denser gas than the brief phase of chemical 
enrichment in gas too diluted to significantly excite the 
HCO+(1-0) line. 

At the opposite, density does not weight the ob- 
servations in absorption. The warmest and richest 
phase of the evolution, may have been detected by 



Liszt fe Lucas (2000)| as the broad and weak absorption 
seen in the direction of B0355+508 over the whole veloc- 
ity extent of the H I absorption line. 

More detailed comparison between abundances in- 
ferred from observations in absorption and in emission in- 
volves time and space averages which cannot be achieved 
at this stage of the study. The present observations have 
been obtained on one single ensemble of such dissipa- 
tive structures, possibly generated simultaneously. Thus, 
space/time average is likely more simple than that re- 
quired to interpret absorption measurements that sample 
depths of gas at parsec scales. 

5.2. Further support for non-equilibrium chemistry: 
"warm" molecules abundance ratios 

The present work stresses the need for a non-equilibrium 
mechanism to explain the large abundances of HCO + de- 
tected in gas of moderate density and low shielding from 
the ambient ISRF. In this scenario, the gas is driven lo- 
cally past a high temperature threshold by an impulsive 
release of turbulent energy. Once the gas has escaped the 
region of large dissipation, it cools down and condenses 
and the chemistry evolves accordingly. 

Our results are independent of the details of 
the processes at the origin of the warm chem- 
istry^ |Flower fc Pineau des Forets (1998) and 

Gredel, Pineau des Forets & Federman (2002) have 
shown that slow MHD shocks of velocities V s ~ 10 km 
1 reproduce the OH/HCO+ abundance ratios observed 



in the diffuse medium ( |Lucas fc Liszt (1996)| ). In MHD 
shocks, as in magnetized coherent vortices, gas is tem- 
porarily heated to temperatures as large as 10 3 K. The 
triggered warm chemistry is characterized by coexisting 
large abundances of OH, H2O, HCO + and CH + , water 
and HCO + being daughters of OH and CH + , respectively 
(JFPF). Their time-dependent evolution in the relax- 
ation phase, is displayed in Fig. ^2 for one of the cases 
discussed in Section 4 (initial density Uh,i = 100 cm" 3 , 
Td = 10~ 21 ergcm _3 s _1 and Ay = 1). The abundance 
ratios along this evolution are interesting in their marked 
differences (Fig. 112(1 . The H2O/OH ratio varies by less 
than a factor 3 along the time-dependent evolution 
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Fig. 11. Time-dependent evolution of (top) the OH (thin 
solid), H 2 (thick solid), HCO+ (dashed), and CH+ (thin 
dashed) abundances, (middle) the kinetic temperature 
and (bottom) the total hydrogen density, along the isobaric 
cooling of a parcel of gas formerly enriched chemically in a 
burst of dissipation of turbulent energy. The initial condi 
tions here are njj,i =100 cm" 3 , 
and the shielding from ISRF Ay=l. 



T d = IQ- 21 ergcm" 3 s" 1 



and up to 10 3 yr after the beginning of the relaxation, 
stays close to the values found in diffuse molecular 
clouds by the SW AS satellite ( |Neufeld et al. (2002)[ 



Plume et al. (2004)| ). The OH/HCCF ratio varies by 
more than a factor 10 along the same evolution. We note 
however that its variations bracket the average value 
OH/HCO+- 30 found by |Lucas fc Liszt (1996)| in their 
random sampling of the diffuse ISM provided by their 
absorption lines survey towards extragalactic sources. In 
contrast, the CH + /HCO + ratio varies very little, before 
10 3 yr only. The value of the ratio, though, critically 
depends on the density (JFPF). This follows from the 
fact that in our model of warm chemistry, HCO + is not 
only a daughter molecule of CH + . It also forms via other 
routes such as CO + + H 2 — > HCO + + H, triggered by 
the endothermic charge exchange between H + and O 
(AE/k = 227 K), and H+ + CO -> HCO+ + H 2 , also 
slightly enhanced in warm chemistry because both CO 
and H;j~ have increased abundances in the warm layer of 
the vortex. This route may even be a dominant formation 
scheme for HCO + , if the large abun dance of Hjj~ detected 
on the line of sight towards £ Per by McCall et al. (2003) 
happens to be the rule in the diffuse ISM. 

There is also an efficient channel to form C 2 H in the 
warm chemistry. It proceeds via the endothermic reaction 
(AE/k = 1760 K) CH + H 2 ^ CH 2 + H, followed by 
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Fig. 12. Time-dependent evolution of abundance ratios, 
OH/HCO+ (thick), H 2 0/OH (thin) and CH+/HCO+ 
(thick dashed), along the same isobaric cooling as in 
Fig. EH 



a series of fast reactions: 



CH 2 



C 2 H+ + H, the 



hydrogenation of C2H + into C2H2 the recombination of 
C2HJ and terminating with photodissociation of C2H 2 
which provides C2H. 

These results are encouraging and suggest that impul- 
sive chemical enrichment of the CNM by turbulence dis- 
sipation bears observable signatures up to several 10 3 yr 
after the end of the enrichment. 



5.3. Large scale signatures of small scale activity. 

Other signatures, independent of chemistry, have been 
found for the existence in the CNM of gas at a 
few hundred Kelvin. They comprise collisional excita- 
tion of H 2 far from UV sources susceptible to gener 
ate fluorescence i Gry et al. (2002)| Lacour et al. (2005) 



Falgarone et al. (2005) 1. In the former cases, H2 is de- 



tected in UV absorption against nearby late B stars, in 
the third case, the pure rotational lines of H 2 have been 
observed along a long line of sight across the Galaxy avoid- 
ing star forming regions. The rotational temperature of 
the warm H2 in these data is 276 K. In all cases the frac- 
tion of warm H2 in the total column of gas sampled is the 
same, of the order of a few percent. It is the same fraction 
as that required in the Solar Neighbourhood to reproduce 
the large observed abundances of CH + within the CNM 
(JFPF). A few percent of warm gas, for which UV pho- 
tons cannot be the sole heating source, is ubiquitous in 
the cold diffuse medium. 

Last, we would like to stress the spatial scales pre- 
sumably involved in the framework of the existing model, 
where individual vortices have transverse hydrogen col- 
umn densities as small as Ah ~ ? N 



sponding to initial values nn,i = 100 cm 3 and li = 15AU 
(JFPF). Wc illustrate the point with the solutions of case 
C: A(HCO+) « 8 x 10 n cm- 2 , A(HCO+) « 2 x 10~ 9 , 
tih ~ 500 cm~ 3 . The total hydrogen column density of ob- 
served material is thus A H = A(HCO+)/A(HCO+) = 4x 
10 20 cm -2 , corresponding to 2xl0 4 vortices in the IRAM 
beam (here we assume that the isobaric evolution occurs 
at constant mass, so that the transverse column density of 
vortices remains the same). For the gas at 500 cm -3 , the 
transverse size has a thickness of about 3 AU. The ensem- 
ble of the 2xl0 4 vortices, probably braided together, as 
vortex tubes do e.g Porter, Pouquet &; Woodward (1994)| 
would form an elongated structure about 450 AU thick, if 
the vortices are in contact with one another. The IRAM 
beam at the distance of the sources is a few times larger 
and the size of the structure seen with the PdBI is of this 
order of magnitude (Paper III). This rough estimate just 
shows that the existence of chemical inhomogeneities at 
AU scales is not ruled out in the diffuse ISM. It would 
certainly help understanding the remarkable similarity of 
the OH and HCO + absorption line profiles reported by 



Liszt fc Lucas (2000) 



10 cm corre- 



Chemical and thermal relaxation timescales are so 
short compared to the dynamical lifetime of molecular 
clouds (of at least a million year) that all the stages of 
the above evolution should co-exist along a random line of 
sight. At present, the receiver sensitivity and the telescope 
beam sizes provide, when observed in emission, a signal 
which is the space and time average of myriads of small 
unresolved structures, likely caught at different epoch in 
their evolution. Molecular line observations in absorption 
in the UV, visible and now submm domains are much more 
sensitive and have presumably already started to sample 
this apparently untractable complexity. 

6. Summary 

This work is part of a broad study dedicated to those 
singular regions in the environment of low mass dense 
cores, characterized by non-Gaussian velocity shears in 
the statistics of the turbulent velocity field of the clouds. 
These regions have been proposed to be the sites of inter- 
mittent dissipation of turbulence in molecular clouds, and 
the large local release of non-thermal energy in the gas 
there has been shown to be able to trigger a specific warm 
chemistry, in particular to explain the large abundances 
of CH + , H2O and HCO + observed in diffuse gas. 

We report here the first detection of HCO + (1-0) line 
emission focussed on two of these singular regions. The 
HCO + abundances inferred from the data cannot be ex- 
plained by steady-state chemistry driven by UV photons: 
the gas temperature is too low and the HCO + abundances 
formed are too small by more than an order of magni- 
tude. As a follow-up of the model proposed by JFPF, we 
have computed the time-dependent evolution of the spe- 
cific chemical signatures imprinted to the gas by the short 
bursts of turbulence dissipation (a few 100 years), once 
the gas has escaped the layers of active warm chemistry. 
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We show that the signature of the warm chemistry sur- 
vives in the gas for more than 10 3 yr and that the den- 
sity (a few 100 cm -3 ), the temperature (a few hundred 
Kelvin) and the large HCO + abundances (2 x 10~ 10 to 
10~ 8 ) inferred from the observations, can be understood 
in the framework of this thermal and chemical relaxation, 
while steady-state chemistry driven by UV photons fails 
by more than one order of magnitude. 

The large abundances detected in the present work, 
and probably the many signatures of a warm chemistry 
now found in diffuse gas, are thus likely due to the slow 
chemical and thermal relaxation of the gas following the 
large impulsive releases of energy at small scale due to 
intermittent dissipation of turbulence. The detailed con- 
frontation of our model with available data relevant to 
the existence of a warm chemistry ongoing in cold diffuse 
gas involves complex time and space averaging, but abun- 
dance ratios among H2O, OH, HCO + and CH + already 
suggest that such non-equilibrium processes at AU scale 
are at work in diffuse clouds. 
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